On Statistics and 1/f Noise of Brownian Motion 
in Boltzmann-Grad Gas and Finite Gas on Torus. I. Infinite Gas 
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An attempt is made to compare statistical properties of self-diffusion of particles constituting 
gases in infinite volume and on torus. In this first part, equations are derived which represent 
roughened but solvable variant of the collisional approximation to exact BBGKY equations. With 
their help, statistics of Brownian motion in infinite gas is considered, under the Boltzmann-Grad 
limit, and shown to be essentially non-Gaussian, involving I/f fluctuations in diffusivity. 
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I. INTRODUCTION 

In p| (see also Part VII and then in 0]) it was argued 
that kinetics of statistically nonuniform (though may 
be thermodynamically equilibrium) gas does not reduces 
to the Boltzmann equation, even under the Boltzmann- 
Grad limit ( v — > oo , S — > , v5 2 = const , with v be- 
ing concentration of gas particles and S radius of their 
repulsive interaction). Moreover, just this limit makes 
clear why the Boltzmann's "Stosshalansatz" ("molecu- 
lar chaos" hypothesis) fails and the Boltzmann equation 
becomes inadequate when the one-particle distribution 
function F%(t, r, v) really depends on spatial coordinate 
r . Instead, an infinite chain of kinetic equations (re- 
produced in Appendix B) arises, which involves, in ad- 
dition to F\ , specific n-particle distribution functions 
F n (t,r,Vi,...,v n ) representing local (probability) densi- 
ties of n-particle chains of collisions (strictly speaking, 
impact parameters and/or other internal parameters of 
the collisions also are arguments of the F n ). The resulting 
long-scale statistics of self-diffusion (Brownian motion) 
of gas particles displays low-frequency fluctuations (1/f- 
noise) of diffusivity, which correspond to not the Wiener 
process but essentially different, non-Gaussian and non- 
Marcovian, random walk (thus confirming earlier phe- 
nomenological considerations 010,0)- 

In contrary to that conjectures, the mathematical the- 
ory of dynamical systems with "hyperbolic behavior" (ex- 
ponential instability and mixing of phase trajectories), 
in particular, such as N hard balls in a container or 
on torus (most interesting for me), succeeds in proof of 
various ergodic, Bernoullian and Marcovian properties of 
these systems (see e.g. 00,13 and references therein). 
In respect to hard balls on torus, this implies, undoubt- 
edly, that their Brownian motion along it reduces (at suf- 
ficiently rough time scales) to the usual Wiener process 
(one can easy make sure of this by numeric simulation of 
3 hard balls on 2-D torus). 

It would be interesting to reconcile these rigorous re- 
sults with above mentioned non-rigorous results on the 



Boltzmann-Grad gas (which can be placed on torus as 
well). There is no contradiction, in principle, if maximal 
time scale, r ng , which still perceives "anomalous" non- 
Gaussian and non-Marcovian features of Brownian paths 
in iV-particle system (or minimum time scale which still 
suppresses and loses these features) , is a growing function 
of N , for instance, 
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where Tf(N) is mean free flight time (we take in mind 
that 1 <C N < oo , of course). 

I can suggest heuristic arguments in favour of this hy- 
pothesis. Indeed, consider Brownian path of some par- 
ticle during some time interval At . The corresponding 
fragment of the path consists of ~ n = At/rf collisions 
and free flights and is fully characterized by Na ~ 22?n 
details (2? denotes torus dimension), namely, velocity 
vectors Vj and flight times Tj ( j = 1 -j- n ) and be- 
sides (2?— l)-dimensional impact parameters of collisions. 
Compare this quantity with a number N c of initial con- 
ditions (parameters of the gas state at the beginning of 
the interval) which might have an influence on this frag- 
ment. It is natural to expect that statistics of the Brow- 
nian displacement AR S = V\T\ + ... + v s t s ( s = 1 -j- n ) 
is sensitive to the ratio Nd/N c . If Nd/N c -C 1 then 
a particular realization of AR S is hardly able to try all 
principal possibilities allowed by initial conditions, there- 
fore, it can be essentially non-ergodic, in one or another 
sense. In opposite, if Nd/N c ^> 1 then any of qualita- 
tively different possibilities will be sampled many times, 
which gives grounds for good ergodicity. 

To estimate N c , notice that outside particle coming to 
j-th collision ( j = 1 n ) previously had ~ j collisions 
with other particles, which in their turn also had colli- 
sions before, and so on. Therefore total number, N p , of 
particles whose initial states influence our fragment can 
be estimated as N p ~ (n + l) n+1 /(n + 1)! - e"/V™- 
At At/rt — n > In N , in fact, initial state of the 
whole gas is important, N p — N . Correspondingly, 
N c = 2VN p - 22? min(AT, e n / y/n) . 

In case of the Boltzmann-Grad gas N — oo , and 
Nd/N c ~ n 3//2 e _n is certainly very small. Hence, non- 
ergodicity of AR n (more precisely, "anomalous" behav- 



2 



ior of fourth and higher-order cumulants of AR n 0,0) 
is not surprising. At finite TV (and not too short time in- 
tervals, At > T f In N), N d /N c = M/NTf . Hence, the 
border between ergodic and non-ergodic behavior lies at 
At/r/ ~ N . That is nothing but our hypothesis Q). 

In view of above reasonings, a system of large fi- 
nite number of particles on torus is of special interest. 
Unfortunately, even in the much simpler case of infi- 
nite Boltzmann-Grad gas the chain of kinetic equations 
discussed in 0, 0] seems not solvable somehow except 
their more or less caricature roughening. Therefore, in 
the present paper I confine myself by two unpretentious 
tasks. In this first part of the paper, the mentioned equa- 
tions are subjected to even more grotesque simplification 
but allowing for their formally exact solution. At that, 
their principal structure is preserved, therefore the so- 
lution is must be quite meaningful. In the next second 
part, modification of the roughened equations for finite 
N will be derived and adapted to torus. Their analysis 
indeed demonstrates existence of characteristic time r ng 
which both formally and essentially corresponds to (IJ. 
Other results will be commented in final Conclusion. 



II. EQUATIONS OF SELF-DIFFUSION IN 
INFINITE GAS 



Our consideration will be based on equations of the 
collisional approximation (see Appendix B and @|), 
first of all, on equations i|B13(l adapted for thermody- 
namically equilibrium but statistically non-uniform gas 
when Boltzmannian collision operator is replaced by the 
Boltzmann-Lorentz collision operator @ . 

Recall that equations l|B13(l deal with n-particle dis- 
tribution functions F n (t, R, V\...v n ) which describe close 
encounters of n particles (with pair collisions between 
some of them) clustered in vicinity of given point R . The 
closeness means that the particles belong to a "collision 
box" defined by \r$ — R\ < d, where d is suitable scale 
much greater than interaction radius 5 but mush smaller 
than dimension ~ (n/i>) 1 /' D of volume typically needed 
to contain n particles ( v and D are mean number of 
particles per unit volume and space dimension of gas, re- 
spectively). At that, the general position n-particle dis- 
tribution functions (DF) F n (t, r\...r n , Vi...v n ) , for par- 
ticles out of any encounter, are factored into product 
nj=i F n (t, Tj, Vj) , at least, under the Boltzmann-Grad 
limit v — > oo , 5 — > , A ^ l/vS D ^ 1 = const , where A is 
mean free path. These DF are normalized to gas volume 
£1 : J F n (t, r\...r n . V\. ..v n ) dv dr — f2 n . Below it will be 
convenient to name F n (t, R, v\...v n ) cluster distribution 
functions (CDF) while F n (t,ri...r n ,Vi...v n ) volume dis- 
tribution functions (VDF). 

The stationary solution to the equations l|B13(l is 



F n (R,vi...v r , 



3=1 



F n eq (v) 



where Fq(v) = (27rtig) _ ' D / 2 exp (— v 2 /2vq) is equilibrium 
velocity distribution, — T/M , and T is gas tempera- 
ture. Such DF represent statistical ensemble which tells 
nothing about positions of particles. 

Now, consider more interesting ensemble which tells 
us that one of particles is localized, say, nearby coor- 
dinate origin, with coordinate probability distribution 
W(r) (JW(r)dr = 1, JrW(r)dr = 0). To write 
this information on VDF, at first suppose that gas vol- 
ume and total number of particles N = vU, arc finite. 
Then the highest VDF can be written as 



N 



F N = F°«(v)^'£w(r j ) 

3=1 

According to relation F n — J F n+ \ dv n+ idr n+ i/Q , the 
corresponding lower VDF are 



F n = F«(v) I 1 + - 



(2) 



Tending £1 to infinity and choosing the distribution 
W(r) to have a suitable width (e.g. ~ A), it is easy to 
see that corresponding CDF can be approximately rep- 
resented by 



F n (R, Vl ...v n ) = F„ e » 



1 



- W(R) 
v 



(3) 



In the Boltzmann-Grad limit (when size of any "collision 
box" turns into zero in units of A ) this is formally exact 
expression. 

Hence, in order to investigate self-diffusion of initially 
localized particle, we should solve equations l|B13p start- 
ing from statistically (but not thermodynamically) non- 
equilibrium initial conditions presented by @. Obvi- 
ously, contribution of the first term on right-hand side of 
©, that is uniform background of densities of particles 
and clusters, will be timc-indcpcndcnt. Therefore we will 
omit it from solution, together with constant multiplier 
\Jv from its time-dependent part, but keep designation 
F n for the rest of it. Thus we come to the problem 



dt dR 



n „ 

Aj / F n+ i dv n+1 

3=1 J 



(4) 



for infinite set of functions F n (t, R, v\...v n ) under initial 
conditions 

n 

W(R)Y[F (vj) (5) 

3=1 

with V^ n > = (vi + ... + v n )/n and Aj being Boltzmann- 
Lorentz operator which acts onto j-th velocity. 

Eventually we are interested in spatial distributions 



W n (t, R) 



F n (t, R, v 1 ...v n ) dv 1 ...dv n 
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In view of their definition and normalization of VDF, 
their statistical sense is quite clear: W n (t, R)/N is den- 
sity of probability that at time moment t > an n- 
particle cluster centered at point R contains (or at 
n = 1 coincides with) those unknown particle which at 
t = was localized near point R = . Naturally, the 
larger is cluster the greater is this probability. Equiva- 
lently, W\ (t, R) is conditional density of probability that 
well-known particle, at t = positioned near R = , 
later will be registered at given point R . In other words, 
Wi(t,R) presents probability distribution of Brownian 
displacement of any gas particle during time interval t . 

Of course, clusters (encounters) of n > 1 particles are 
ephemeral objects whose lifetime ~ d/vo is much smaller 
than mean free flight time r/ . But they form constantly 
living "gas of encounters" under description by second 
and higher of equations . The latter have no principal 
difference from the first equation since it also describes 
statistics of gas but not individual particles. 



III. ROUGHENED DESCRIPTION OF 
BROWNIAN MOTION 

1. In the Boltzmannian kinetics, instead of Q one 
has the single equation for one-particle DF F\(t,r,v) , 



dF\ 
~dt 



dF\ 
dr 



vKFx 



(6) 



where the Boltzmann-Lorents operator A tends ve- 
locity distribution of gas particles to equilibrium one, 
Fq(v). As above, we can treat F\{t,r,v) as (small) non- 
stationary non-uniform addition to uniform background 
due to a priori information about initial position of one 
of particles. After solving © with corresponding initial 
condition F\(t — Q,r,v) = W(r)Fo(v) , the distribution 
W\{t,r) — jFi(t,r,v)dv will present us information 
about a law of Brownian motion of the selected particle. 

Undoubtedly, in respect to large time scales, S>t/, 
that will be nothing but the Gaussian law, which can be 
expressed in several equivalent forms: 



lim Wi(s 2 t : sr) = {ATiDt)-' D/2 exp {~r 2 /ADt) , 

s — >oc 

lim Ex(t/s, y/sk) — exp (— Dtk 2 ) , 

s— >o 

lim s$i(sp , Vsk) = l/(p + Dk 2 ) 



(J) 

(8) 
(9) 



Here Si(t, fe) is characteristic function of Brownian mo- 
tion and Si (p, k) its Laplace transform: 

3i(t, k) = J exp(ik-r)Wi(t,r)dr , (10) 

*i(p, fe) = / exp(-pt)Z 1 {t,k)dt (11) 
Jo 

Such long-time asymptotic is insensitive to concrete 
form of operator A . Therefore, let us learn to obtain it 



in most simple way, and then apply this our experience 
to equations 

The Boltzmann-Lorentz equation directly yields 



0Wi(t,r) dJ^r) 



dt 



dr 



(12) 



J x (t,r) = J vF 1 {t,r,v)dv = v(t,r)W 1 (t,r) , (13) 

d{Jl{t dt r))a =-Jj-j 'vv a F 1 dv + v J v a AF 1 dv (14) 

Here, evidently, J\(t,r) represents probability flow and 
v(t,r) conditional mean value of velocity under given 
position of the particle. 

Notice that soon after start there is some correlation 
between current velocity of the particle and its displace- 
ment from initial position. But then, at t 3> Tf , this 
correlation decreases down to zero for any thermodynam- 
ically equilibrium, and hence symmetrical, random walk 
(otherwise it would behave like ballistic flight). In other 
words, the conditional mean velocity is much smaller 
than above defined mean-square thermal velocity vq : 



v(t,r)/v -> 



(15) 



for any possible realization of the walk (more precisely, 
one could write v(t, r)Wi(t, r)/v W\{t, 0) — * ). Corre- 
spondingly, Fiit^r^v) factorizes into Wi(t,r)F (v) at 
t ^> Tf , and integral in the first term of l|14|) is dominated 
just by this factorized asymptotic: 



d_ 

dr 



f d 
/ vv a Fx(t,r,v)dv -> v 2 — Wi{t,r) 
J or a 



(16) 



In opposite, second integral in (|14fl . as well as the 
probability flow Ji(t, r) , is determined just by the weak 
velocity-displacement correlation. Transform it as 

v J v a AFi(t, r, v) dv = v J Fx(t, r, v) [A^ v a ] dv , 

where symbol ' stands for transposition. In view of time- 
reversal symmetry of collisions, v a is an odd function 
of v . Therefore the latter integral is an odd function 
of the conditional mean velocity v(t, r) . At f > t/ , 
in view of (|15f) . is turns into linear function of v{t,r) 
and thus of J\{t, r) . The coefficient follows from the 
asymptotic factorization of F\ , and the result is 



v J v a AF 1 (t,r,v)dv = -j(Jx(t,r)) a , (17) 

Thus, velocity of particle is excluded from our consid- 
eration. Combining ifT^jl. ifTljl. with (fTTjl . and using 
V for spatial derivative, we come to approximate equa- 
tions which deal with particle displacement only: 



dW 1 



r\ T 

-V-Ji , = -«gVWi-7Ji (18) 
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Corresponding Laplace transform lllj l of the character- 
istic function I|1U|) looks as 



*i(p, fc) = Si(0, fc) 



P + 7 



(19) 



and evidently possesses Gaussian asymptotic (0, with 
diffusivity D = Vq/'j. At that, velocity relaxation rate 
7 can be identified with tJ 1 . 

Let us compare (|19fl with formally exact answer when 
A is Fokker-Planck operator 0,0 governing random pro- 
cess with Gaussian velocity fluctuations (the Ornstein- 
Uhlcnbeck process), 



gifo fc) 
Hi(0, fc) 



v 2 k 2 



p + 7 + 



2i> 2 fc 2 



p + 2 7 + 



This is infinite continued fraction with tivq fc in nom- 
inators and n/j in denominators (n = 1, 2, ... ). As for 
equations (|18|l and (|19|) . it can be shown that formally 
they correspond to a random process with dichotomic ve- 
locity fluctuations. However, from the point of view of 
large time scales both random walks are identical. Re- 
alistic A's account for multi-state velocity hopping but 
lead to same random walks. This example gives hope 
of that similar approach to system @ will imitate its 
long-time behavior without principal losses. 

2. In addition to distributions W n (t,R) , introduced 
in Sec. 2, consider the "flows" 

J n {t,R) = J V {n) F n (t : R,v 1 ...v n )dv 1 ...dv n 

In essence, equations (@J describe the same physics as 
equation Hence, all above reasonings about decay (at 
t Tf ) of correlations between velocities and positions 
of particles remain valid. And all they can be confirmed 
by direct analysis of @. Therefore, by analogy with (|16|) 
we can write 



V^V^F n (t 7 R, Vl ...v n )d n v 
vWyWf^^) d n v V W n (t, R) 



= —V a W n (t,r) (20) 
n 

(here and below d n v = dv\...dv n and V replaces 
d/dR). Analogue of formula (fT7|l is 



v / F n+1 



« n 

3=1 
1 " 

n ^ J ■ 

j=i 

, n+l 



v \ F n+1 



d n+1 v 



3=1 



d n+1 v 



fJ n+1 (t,R) (21) 



with the same relaxation rate 7 as above. We took into 
account that Ajl = , therefore Aj v m — at j ' ^ m , 
and that, naturally, all F n (t, R, v\... v n ) can be treated 
as fully symmetric functions of velocities. Besides, now 
each of F n involves its own conditional mean velocity 
v n (t,R) . This derivation becomes quite transparent in 
the case when v is eigenfunction of A : Av = — jv . 
Consequently, we came to equations 



dW n 
dt 



dJ n 

dt 



= -V • J n , 
2 



^VW n - 1 J n+l 



(22) 
(23) 



with initial conditions following from JSJ): 

W n (t = 0,R) = nW(R) , J n (t = 0,R) =0 (24) 

At that, since a shape of W(R) will be of no importance, 
we can simply put on W(R) — S(R) . 



IV. STATISTICS OF SELF-DIFFUSION AND 
LOW-FREQUENCY FLUCTUATIONS IN 
DIFFUSIVITY 

1. Applying to l|22l) - (|23|l Fourier and Laplace trans- 
forms, one can easy obtain formal solution to the problem 
(l22H - l|24f) in the form of infinite series: 



, 2 00 f 



s—n m—n 



7 



p m + 7_Dfc 2 jp 1 



(25) 



Here <fr ra = $n(p, fc) are defined analogously to (JUJ, and 
D = ^0/7 ■ I11 particular, at n = 1 we have 



p$i = 1 



P 



) 






p. 



(26) 



with B(x, y) being beta-function. After that, using 
standard integral representation of this function, we can 
transform the latter series to compact integral form: 



p$i = 1 - Z 



{i + jx/ P y 



dx , Z^ 1 -^- (27) 



Let us discuss statistical properties of Brownian motion 
corresponding to this formal result. 

2. First, consider what is long-scale asymptotic of our 
random walk, if interpret its asymptotic similarly to Q- 
@. In the limit like in §5J, equation yields 



p lira s$i(sp, y/sk) = 

s—>0 



(28) 



1 - 



Dk 2 



1 - 



P Jo 
Dk 2 



exp 



Dk 2 



dx 



(Dk 2 ) 2 



P 



(l + x) 2 

- X 

dx 



Dk 2 +px 
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The last row is obtained with the help of integration by 
parts. Integration by parts "in opposite direction" gives 

lim s^(sp, sTsk) = 2 r eX P ( ~ Jfc ' T) dr (29) 
s^o J Q (l+pry 

Performing here inverse Laplace transform, we obtain 



lim Ei(t/s, \fsk) — 

s^O 

/ exp(-xDfc i) exp — — 
Jo \ x) x A 



(30) 



with Hi defined by H10JI . Then the inverse Fourier trans- 
form of i|3U|) produces 



lim W 1 (s 2 t,sR) ■■ 

s — >oo 

T(2 + V/2) 



(31) 



(4irDt) v / 2 (1 + R 2 /ADt) 2 + v / 2 



This probability distribution describe random walk 
which has the same characteristic scaling property R 2 ~ 
t as Gaussian random walk (Wiener process), but at 
that obeys extremely non-Gaussian statistics. Formula 
(|3U[) shows that this walk can be represented as Wiener 
process with random but time-independent diffusivity D 
distributed with probability density w(D/D)/D, where 



w{x) 



' exp (— l/x) . That is absolutely non-ergodic 



process. Besides, as it is evident from (|31|) . the fourth 
and higher statistical moments of this process are infi- 
nite. Apparently, non-Gaussian and non-ergodic pecu- 
liarities of our random walk was too exaggerated by the 
considered limiting procedure, and we have to return to 
equation l|27l) which demonstrates that in fact all the mo- 
ments are finite. 

3. Let R(t) denotes projection of 2?-dimcnsional 
Brownian displacement R(t) during time t onto some 
fixed direction, and vector k is oriented in parallel to 
it. Since the direction is fixed, we can treat k as scalar. 
Then compare expansions of both sides of 1)2 7|l into series 
over k 2 : 



p ~y ( 2n r- Jo 



,-pt lr>2n 



p ±i 



p 

jDk 2 

9 

P 



1 



(R 2n {t))dt = 



[- ln(l -x)] n 

(n-i)i y (i + jx/ P y 



■ dx 



At p/j <C 1 , which correspond to large time, the com- 
parison yields 



9 D 

e~P t (R 2 (t))dt = -r , (32) 
o P 

e-r* (R 4 {t))dt = Infl + ^ , (33) 

p 3 V p) 



-pt ( R 2n^ dt = 



(2n)\ C(n 



D 2 v 2r 
I) V ° 



(n > 2) (34) 



Here ((■■■) is the Riemann zeta function, and we used 
the integrals 



1 [-ln(l-x)] n , ... , . 

l - ^ '-^-dx = nl((n) (n>l), 

x z 

f 1 ln(l-x) , hx(l + a) 

I dx —— — 

J (l + ax) 2 ' a(l + a) 

Correspondingly to l|3^ )l -l|3l )l . 

(R 2 {t)) = 2Dt , 
(R 4 {t)} = \2D 2 t 2 In 7* , 



(35) 
(36) 



{R 2n (t)) = 2n(2n - 1) Q{n - 1) x 



2 /„, j.\2n- 4 



x {my (v t) 



(n > 2) (37) 



at jt 3> 1 . Thus, except mean square of the displace- 
ment, all its higher-order moments are anomalously large 
from the point of view of Gaussian random walk (which 
has (R 2n (t)) = (2n)!! (2Dt) n ). 

Comparing formulas (|35|) - (|37|l and distribution l|31|) at 
T> = 1 , which also describes the projection R(t) , we 
notice that these asymptotic formulas can be agreed with 
probability distribution l|31|) if the latter is supplemented 
by its cut off at \R(t) \ ~ vot . Hence, we can improve l|31|) 
if replace it, in asymptotical sense, by 



Wi(t,R) 



r(5/2) G(R/v t) 
(AirDt) 1 / 2 {I + R 2 /ADtf/ 2 



(38) 



where G(0) = 1 and G(x) is extremely fast decaying 
function at x — * ±oo (this statement certainly follows 
from very slow growth of coefficients on right-hand side 
of (JS3). Thus we come to conclusions as follow. 

(i) In respect to probabilistic characteristics of out ran- 
dom walk, formula l|38|l practically coincides with Q3ip. 
that is (|31|l is quite satisfactory. However, in respect to 
higher statistical moments, (R 2n {t)) with n > 1 , there 
is not only strong quantitative difference between 131|) 
and H38f) but also at least two principal differences. 

(ii) The first difference is presented by fourth (R 4 (t)) 
moment. Such behavior of (R 4 (t)) as prescribed by 
(|36|l means that our Brownian motion indeed can be 
represented as Gaussian one with random diffusivity. 
However, diffusivity is not a "static" random value (as 
(|31|l would imply) but (as l|38f) implies) a random pro- 
cess with characteristic low-frequency power spectrum 
~ (In/)' 3 // (where / <C 7 is frequency and (3 can be 
both positive and negative). In other words, the diffu- 
sivity fluctuates like 1/f-noise. This theme was already 
discussed not a once in [lj, |2j, |3j, |4j, |fj . 

(iii) The second peculiarity of our Brownian motion 
is presented by moments (R 2n (t)) with n > 2 which 
turn out to be dependent on not only (average) diffusiv- 
ity D but also on the thermal velocity vq , i.e. charac- 
teristics of separate free flights. This formal result can 
be interpreted as evidence of that probability distribu- 
tion of free flight time has a power-law tail. Comparing 
(|36|l and moments corresponding to purely free flight, 
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(R 2n (t)) = {2n)\\ {v t) 2n (on (R 2n {t)) = (v t) 2n for 
dichotomic velocity), we can suppose that probability of 
anomalously long flights comparable with t behaves as 
~ (Dt) 2 /(vot) 4 = l/(7t) 2 (i.e. is inversely proportional 
to square of mean number of collisions during the ob- 
servation time). But, strictly speaking, this issue needs 
in investigation of many-time statistics of our Brownian 
motion. In any case, on the whole it appears to be rather 
rather non-ergodic process. 



RESUME 



dF 2 (9) 

"aT = L F2 
at 



(W 



>L {L > 



dF n 
dt 



= L (n) F 2 + (N-n) 



L 23 )F 3 



(A3) 



(A4) 



In 



where J s ... = J ...dv s dr s and F n = J n+1 F n+1 
the limit of infinitely big container, N — » 00, I — ► 00, 
v = N/il = const (51 = l v ), one must consider non- 
normalized distribution functions Sl n F n . If designate 
them again as F n then (|A4() turns into 



In this paper, previously constructed equations of the 
collisional approximation to kinetics of statistically non- 
uniform gas were applied to analysis of self-diffusion 
(Brownian motion) of particles of infinite-volume gas in 
its thermodynamical equilibrium, under the Boltzmann- 
Grad limit. With this purpose, a roughened version of 
the equations was derived which allowed for their almost 
complete solution. The latter demonstrated that statis- 
tics of the Brownian motion is essentially non-Gaussian, 
not obeying the central limit theorem, and hence essen- 
tially non-ergodic. It looks as if diffusivity of the Brown- 
ian path had no definite value but instead underwent low- 
frequency fluctuations with spectral properties similar to 
that of 1/f- noise and with probability distribution pos- 
sessing power-law long tail. These results are in abrupt 
contradiction to what follows from standard kinetic ap- 
proximations, but agree with Krylov's doubts [Tol ] about 
obligatory ergodicity of mixing dynamics. 



APPENDIX A: BASIC EQUATIONS 

Our start point are the Liouville equation for N par- 
ticles in container < r a < I (a = 1 -r T>): 



OF, 



N 



-(AT) 



dt 



N 

3=1 



(Al) 



9r, 



Fj m , 

l<j<m<N 



L, 



-Vf/fo - r m ) 



d 



d 



dpj dpr, 



with U(r) being interaction potential, pj = Mvj mo- 
mentums, and Fn = F^(t, rx...ri\r, v i...vi\r) normalized 
probability density of phase point of the system, and the 
corresponding Bogolyubov-Born-Green-Kirkwood-Yvon 
(BBGKY) equations 

-^ = - Vl .^- + [ N-l)j^F^ (A2) 



dF n 
dt 



L^Fn 



iF, 



n+l 



(A5) 



The equations l|Al|l - ljA4|l hold also for particles on 
torus with equal dimensions, < r a < I , if only supple- 
ment them with the periodic boundary conditions and 
replace U (r) by U (p) , where p = p(r) is factual dis- 
tance between two particles on torus: p a = r a when 
\r a \ < 1/2 and p a = — sign(r Q )(^— \r a \) when \r a \ > 1/2. 
However, if we want to consider unbounded Brownian 
motion of particles along the torus, we should introduce 
more rich distribution functions which keep count of their 
rotations. This is subject of second part of the paper. 



APPENDIX B: COLLISIONAL 
APPROXIMATION 

Assuming that our gas is sufficiently rarefied, that is 
S/a <C 1, with a — l/N 1 ^ being typical distance be- 
tween neighbour particles and S interaction radius, we 
want to describe it approximately in terms of Boltzn- 
mannian pair collision operators ("collision integrals"). 
With this purpose, basing on the exact BBGKY equa- 
tions (see Appendix A) consider F2(t,r\,r2,vt,V2) at 
1 1*2 — n| < d where {«d«a. Any such configura- 
tion can be interpreted as a snapshot of certain process 
of collision of two particles or, may be, their close en- 
counter without significant interaction. Therefore let us 
name the region \r 2 — r±\ < d 2 "collision bo^ and char- 
acterize states inside it by inner time O of corresponding 
collision plus relative velocity and (T> — l)-dimensional 
impact parameter at its incoming stage. We can choose 
= at perigee of the collision or encounter. If using 
these new variables, (IA3I) can be rewritten in the form 



+ (N-2) (L 13 + L 23 )F 3 



dF 2 = (2) 8F 2 dF 2 

dt ' OR dG 

(Bl) 

with R = (n +r 2 )/2 and V {2) = {v x +v 2 )/2 being 
position and velocity of center of gravity of the pair. At 
that, by definition of the inner time, 



dF 2 
d® 



dF 2 
dri 2 



L12F1 



(B2) 
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where v\ 2 = r 2 — T\ and V\ 2 — v 2 — v\ . 

Any dependence of F 2 on 6 , at least at |0| ~ typi- 
cal duration of collisions r c , does mean that at statistical 
ensemble under consideration a number (or probability of 
realization) of pairs coming into ( O < ) a given sort of 
collisions differs from a number (or probability of realiza- 
tion) of pairs leaving ( O > ) the same collision. Gen- 
erally such is the case, and such difference obstructs col- 
lisional approximation, because collision operator must 
conserve probability (just as collision conserves number 
of particles). Or alternatively we should neglect this dif- 
ference and postulate that 



dF 2 
90 



when r 12 S collision box 



(B3) 



111 CBTJ, thus ensuring conservation of probability in the 
course of collisions. The destination of l|B3[l is correct 
transition from separate dynamic states to coherent col- 
lision as the whole. 

According to l|B2|) . relative displacement of colliding 
particles is included into collision. Consequently, it drops 
out of outward characterization of collision as the whole, 
and the latter obeys the equation 

BF r)F P 

~dT = _y(2) ' ~di + {N ~ 2) J 3 {Ll3 + L23)F3 (B4) 



At that, equation (|B2|I . which is probabilistic equivalent 
of Hamilton dynamics inside collision box, 



dr 12 dv 12 
v 12 , Ai- 



de 



dO 



dU{ ri 
dr 



12 



(B5) 



together with i|B3(l serves for transformation of last term 
in l|A2(l into usual Boltzmann collision integral Q, Q] . 

Here we detect principal difference from Boltzmann's 
theory: equation (|B4|) certainly forbids Boltzmann's 
Stosshalansatz (except trivial situation when F± is ex- 
actly spatially uniform). "Molecular chaos" fails in col- 
lision box, just where it was most wanted. However, 
nothing prevents "molecular chaos " outside collision box, 
and it really holds there (at least in the Boltzmann- Grad 
limit). 

Clearly, the function F 2 governed by equation l|B4|) 
presents probability measure (ensemble-average density) 
of one or another sort of collisions (differentiated e.g. by 
their impact parameters and relative velocity Vi 2 ). No- 
tice that it becomes normalized to unit, in respect to 
I?( 2 ) and both velocities, if multiply F 2 by gas volume 
fi . The first term on right-hand side of (|B4|) describes 
drift of collisions in spatially non-uniform gas while the 
second their redistribution over sorts due to intervention 
of a third outer particle. 

Hence, next we should consider clusters of three close 
particles belonging to a reasonable "collision box", e.g. 
|r*j — Tj\ < d, and again connect instant dynamic states 
inside it into 3-particle collision as a single whole. Of 
course, that is not literally 3-particle collision but sooner 
3-particle encounter with possibility of pair collisions 



among some of the particles. Again, distribution of such 
3-particle events over their sorts can change under in- 
fluence by a fourth exterior particle. Continuing such 
scheme to n-particle encounters, we decompose the n- 
particle Liouville operator similarly to i|Bljl: 



L {n) F n = -V {n) ■ 



dF n 8F n 
dR dQ ' 



(B6) 



with R= (r 1 + ... + r n )/n , V< n > = (v 1 + ... + v n )/n , 
and being inner time of an encounter. Then again it 
is necessary to postulate conservation of probability in 
collision box: 

dF n 

~~qq = when — rj G collision box , (B7) 



(B8) 



which simplifies (|A4|) into equation 



9Fn = _ v (n)9F n 



ot 



dR 



-(N-n) 



/ ^ Lj n+ iF n+ i 
Jn+1 - =1 



similar to l|B4(l . At the same time, (|B7(1 helps to trans- 
form the last term in similar equation for F n _i into sum 
of n — 1 collision integrals. 

Thus we come to a chain of coupled kinetic equations: 



dF„ 
dt 



dFn 
dR 



+ (N -n)Y^S jn+1 F 



3=1 



n+l 



(B9) 



where n = 1...N and Sj„+i is Boltzmannian collision 
operator describing collision between j-th particle from 
the n-particle cluster and one more " [n + l)-th" parti- 
cle from its exterior. Functions F£+i represent F n+ i 
at boundary of (n + l)-particle collision box, more pre- 
cisely, at those part of the boundary which correspond 
to incoming stage of these collisions. At that, differences 
~ d < a between centers of gravity of successive colli- 
sion boxes are neglected because of smallness in compar- 
ison with the mean free path A ~ a(a/5) v ~ 1 . In the 
limit of infinite gas, after normalization to its volume, 
Sl n F n — > F n , equations i|B9|) turn to infinite chain of 
equations from [HQ], 



dF n 
dt 



dF n 
dR 



(BIO) 



Importantly the boundary functions F^^ in (|B9(1 and 
(|B10p in no way are ruled by these equations themselves. 



In essence, F"^ supply initial conditions to equations 
(|B5|) . as well as similar n-particle equations, and must be 
determined independently, by lacing i^+i's with proba- 
bility distributions at exterior of the collision box. But its 
boundary is just border between violation and observance 
of "molecular chaos" (see Italic text above). Therefore, 
here somewhat mediate must be expected. For example, 

F 2 m {t,R, Vll v 2 ) = F{(t,R,v 1 )F{'(t,R,v 2 ) 
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where generally neither F^(t,R,v) nor F"(t,R,v) co- 
incide with the one-particle distribution Fx(t, R, v) . 
Thus Stosshalansatz is valid with respect to velocities of 
colliding particles only but not to their positions ! More 
concrete formulation of such "weakened molecular chaos" 
follows again from conditions of conservation ofprobabil- 
ity (now at the border of "molecular chaos" ) P, Q : 

F$_ 1 (t,R,Vi...v n ,v n+1 ) = Fi(t,R,v n+ i) x 

x J F n+1 (t,R,Vx...v n ,v')dv' , (Bll) 

where Fx is local one-particle velocity distribution, 



Fx(t,R,v) 



Fx(t,R,v) 



J Fx(t,R,v')dv> 



(B12) 



normalized to unit. 

But what about T>(n — 1) — 1 those arguments of F n 
at n > 1 (in particular, impact parameters) which con- 
cretize geometry of n-particle encounters but not written 
out above ? If we made attempt to take them into ac- 
count, our theory would be not simpler than underlying 
BBGKY equations. Therefore we suppose that F n in 
(|B9ll and (|B10I) represent distributions averaged over all 
sorts of the encounters. This roughening of the theory 
looks be caricature, but it can be legalized in many re- 
spects under the Boltzmann-Grad limit ( v — > oo , S — > , 
A ~ l/vS 1 ^^ 1 = const). What is the chief thing, at 
that the theory keeps principal open-chain structure of 



its equations. By this reasons, in Sec. 2 we venture to 
roughen the theory yet stronger. 

Additional simplification appears when we want to in- 
vestigate Brownian motion (self-diffusion) of gas particles 
in absence of their collective hydrodynamic motion, that 
is in thermodynamically equilibrium gas. This purpose 
can be achieved if apply the trick well known in the clas- 
sical kinetic theory of gases , namely, if treat "outside 
(n + l)-th particle" in (IB11I) - (|B12|) as a particle of equi- 
librium thermostat: Fx(t,R,v) — > Fq(v) , where Fq(v) 
is equilibrium velocity distribution (standard Gaussian 
at N — > oo ) . Correspondingly, the Boltzmann collision 
operator Sj n +i transforms into the Boltzmann-Lorentz 
operator Aj , defined by 



S 



.)W (v n+ x)} = \j /'(...r,...i 



jn+X 



Equations ((BTOjl - (JBTTJ take the form [lj 



dFn = _ (n) dK 
dt dR 



A, J F n+1 dv n+1 (B13) 



and similarly equations (|B9|) do change. Thus we find 
ourselves in situation when spatial non-uniformity of dis- 
tribution functions F n has purely statistical origin: it 
describes not a thermodynamic perturbation but infor- 
mation about positions of some particles only (even may 
be a single particle). 
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